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ABSTRACT 

We have conducted B, g, V, and iZ-band imaging in a 45' x 40' field containing part of the high 
Galactic latitude translucent cloud MBM32, and correlated the intensity of diffuse optical light 
S V (X) with that of 100/im emission (100/im). A x 2 minimum analysis is applied to fit a linear 
function to the measured correlation and derive the slope parameter 6(A) = AS 1/ (X)/AS v (100fj,m) 
of the best-fit linear function. Compiling a sample by combining our 6(A) and published ones, 
we show that the 6(A) strength varies from cloud to cloud by a factor of 4. Finding that 6(A) 
decreases as 5,, (100/im) increases in the sample, we suggest that a non-linear correlation including 
a quadratic term of S u (100/im) 2 should be fitted to the measured correlation. The variation of 
optical depth, which is Ay = 0.16 - 2.0 in the sample, can change 6(A) by a factor of 2 - 3. There 
would be some contribution to the large 6(A) variation from the forward-scattering characteristic 
of dust grains which is coupled to the non-isotropic interstellar radiation field (ISRF). Models 
of the scattering of diffuse Galactic light (DGL) underestimate the 6(A) values by a factor of 2. 
This could be reconciled by deficiency in UV photons in the ISRF or by a moderate increase in 
dust albedo. Our 6(A) spectrum favors a contribution from extended red emission (ERE) to the 
diffuse optical light; 6(A) rises from B to V faster than the models, seems to peak around 6000 
A and decreases towards long wavelengths. Such a characteristic is expected from the models in 
which the DGL is combined with ERE. 



Subject headings: diffuse radiation - dust, extinction 
(MBM32) — scattering 



infrared: ISM — ISM: clouds — ISM: individual 



1. Introduction 

A diffuse optical component to the Galactic in- 
terstellar medium (ISM), which is called diffuse 
Galactic light (DGL) or sometimes "optical cir- 
rus," was noticed in late 19 30's. Extensive studies 
(e.g., lElvev fc Roachlll937 



Henvev fc Grccnstcin 



19411 : Ivan de Hulst fc de Jondll96ftlMattUalll979h 
revealed that the DGL is starlight scattered off 



by dust grains in the diffuse ISM which is illu- 
minated by the interstellar radiation field (ISRF). 
The spectrum of the D GL extends t o ultr aviolet 
(UV) wavelengths where Iwitt et all ( 19971 ) made 
extensive studies of the dust scattering properties. 
The DGL distribution is strongly concentrated to- 
ward the Galactic plane. This concentration is 
attributed to the forward-scattering characteris- 
tic of interstellar dust grains coupled with a non- 
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isotropic distribution of stars. 

Over 50 years later, diffuse far-infrared (IR) 
emission, so-called "IR cirrus," was discovered 
i n the IR Astron omical Satellite (IRAS) mission 
(|Low et al.lll984f) . It was quickly shown that the 



IR cirrus is also visible on photographic ma terial 
ride Vries fc Le Pooldll985l:lPalev et al.lll99lh . and 
suggested that the 100/um brightness varies lin- 
early with extinction in a range of u p to ~ 20 MJy 
sr" 1 ( de Vries fc Le Poolei 19851 : Laureiis et al 
1987). Finding good agreement between the struc 
tures seen on the optical images and the distribu- 
tion of the CO emission, IStarkl (|l995h pointed out 
that the optical images can be used to distinguish 
variations in the dust column density from those 
in the molecular column density. These studies in- 
dicate that the IR cirrus is the far-IR counterpart 
of the DGL; far-IR emission in IR cirrus is ther- 
mal emission which follows the dust absorption of 
starlight. It is thus natural to expect a linear cor- 
relation between the DGL and far-IR brightness 
in the optically thin limit. A combination of opti- 
cal and far-IR observations would give a powerful 
tool for investigating the dust properties as we ll as 
the ISRF in diffuse ISM (e.g.,lBeichmarJll987lan d 
references therein; Guhathakurta fc Tyson 1989). 

Since the late 1990s, remarkable progress has 
been made in the optical and far-IR. In the op- 
tical, new wide-field CCD cameras make it pos- 
sible to cover one or more square degrees at a 
time, providing the uniform sensitivity over a 
large area and reducing the systematic uncer- 
tainty by temporal a i rglow variations. In the 
far-IR, ISchlegel et all (|l998h . hereafter SFD98, 
published the IRAS/DUiBE maps, which are 
reprocessed all-sky maps, by combining IRAS 
and DIRBE (Diffuse Infrared Background Ex- 
periment on board the Cosmic Background Ex- 
plorer satellite) data. The IRAS/DIRBE maps 
have DIRBE quality calibration and IRAS res- 
olution and feature removal of point sources as 
well as of artifacts from the IRAS scan pat- 
tern. Recent studies of diffuse optical light benefit 
from the progress using a wide-field CCD camera 
and/or the IRAS/DIRBE maps. Analysis of dif- 
fuse optical light recorded in the Pioneer mission 
against the IRAS/DIR BE lOO^m maps finally led 
Matsuoka et al.l (|201ll ) to report the detection of 
the cosmic optical background (COB). The DGL 
is an important foreground to be subtracted from 



data taken to measure the COB. Thus, studying 
the DGL distribution is indispensable to a better 
understanding of the COB. 

The present paper mainly focuses on the follow- 
ing issues. The first issue is related to variations 
in slope parameters, i.e., AS U (X)/ AS u (100um), 
which is determined by fitting a linear func- 
tion to the measured intensity of diffuse opti- 
cal light S v (A) as a function of 100/im intensity 
S^lOO/im). It has been indicated that slope pa- 
rameters vary from clou d to cloud by a factor 
of 3 - 4. For example, Guhathakurta fc Tvsonl 
(|l989h reported the slope parameters scatter by 
a factor of 3 in their sample of 4 high Galactic 
latitud e clouds, and the slope in a cloud stud- 
ied by IPalev et ajj f)l99lh is 4 times greater than 
that obtained from blank-sky spectra sampled in 
92,000 d irections at high la titude (mostly |6| > 
30°) by IBrandt fc Draind <|2012h . What does 
cause such large variations? Another issue is 
related to extended red emission (ERE) in the 
diffuse optical light. ERE detected in the re- 
flection nebula NGC2327 has a broad emission 
feature with FWHM of 1500A peaking at 7200A 



(jWitt fc Schildl 119881: [Gordon et al.l Il998l) . Car 



rying out a sp ectroscopic survey of high-l a titude 
cirrus clouds, ISzomoru fc Guhathakurta (Il998l) 
found that the cirrus spectrum does not agree with 
that of synthetic scattered E-type galaxy spectra, 
and explained this mismatch by the presence of 
ERE. Their cirrus ERE feature peaked at 6000A 
and is shifted toward short (bluer) wavelength 
as compared with ERE in HII regions (peak ~ 
8000A and reflection and planet ary nebulae (peak 
~ 700 0A). On the other hand, IBrandt fc Draine 
(|2012l ) concluded that their blank-sky spectrum is 
consistent with those expected for their scatter- 
ing DGL models, finding no evidence of ERE in 
the diffuse optical light. For broadband observa- 
tions, some authors favor the presence of ERE in 
the diffuse optical light (e.g., Gordon et al. 1998; 
Witt et alJliopj iMatsuoka et alll201ll). while the 



others suggest the absence of ERE (jZagurv et al 
19991 : IZagurvll200fih 



In the course of our program to constrain the 
COB, we have carried out broadband imaging of a 
45' x 40' field which contains a translucent molec- 
ular cloud at high Galactic latitude. In section 2, 
we describe the observations and the data process- 
ing details, such as flux calibration, masking, and 
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aperture correction. Correlation analysis of the 
intensity of diffuse optical light against the 100/im 
emission is presented in section 3. In section 4, we 
use our measurements along with similar measure- 
ments in the literature to discuss issues affecting 
the correlation, and consider implications for dust 
scattering DGL models. A summary appears in 
section 5. 

2. Observations and Reduction 
2.1. Observations 

The target field was selected to have the follow- 
ing properties: (1) it is optically thin or translu- 
cent at optical wavelengths, and is characterized 
by a sufficient contrast at 100/^m to be distin- 
guished from the 100/^m background; (2) it is not 
crowded with stars and does not include bright 
stars. The first criterion is needed in order to 
see clearly the correlation between optical 100/im 
emission. The second is needed in order to mini- 
mize the contamination from stars, especially from 
the extended halos of bright stars. Using these 
criteria, we have selected the southern part of 
MBM32, covering - 1° x 1?5. This is a dif- 
fuse, translucent molecular cloud with extin ction 
ofA v = 0.3-0.5 (|Magnani fe de Vrieslll986l). and 
locat ed at a distance of 120 pc ( Magnani et al.l 
1996h at Z=147?2, 6=40?7. The surrounding sky 
is not a pure blank-sky; thermal far-IR emission 
from interstellar dust was d etected. Ba s ed on the 
IRAS/BIRBE 100/im map. IWitt et al.l (j2008l) es- 
timated Ay — 0.6 for the cloud and Ay = 0.25 for 
the surrounding sky. The 100/im intensity in the 
field ranges from 1.4 to 7.0 MJy sr -1 . 

Observations were made at Kiso observatorj0 
on dark nights in 2011 February, April, and 2012 
February. The Sun and the Moon were more than 
35 degrees below the horizon when the data were 
acquired. The optical data were acquired using the 
105cm Schmidt telescope with t he 2KCCD camer a 
which has 2048 x 2048 pixels (|ltoh et al.ll200ll) . 



The field of view of the camera is 50' x 50' with a 
pixel scale of 1'.'5 pixel -1 . The average seeing was 
~ 3'.'3. The area of high-quality is approximately 
45' x 40', smaller than the field of view of the 
camera due to using dithering. 



In Figure 1, the observed 45' x 40' field is in- 
dicated with a white square superimposed on the 
IRAS /DIRBE 100/im map. Images were obtained 
at four photometric bands, B, g, V, and R. In order 
to reduce the effects of large-scale non-uniformity 
in the flat-fielding, we observed the field at two dif- 
ferent telescope attitudes by flipping the telescope 
around the right ascension axis in such a way that 
the difference in hour angle is 180 degrees, thereby 
rotating the image plane on the detector by 180 
degrees. In each telescope attitude, the field was 
observed three times per photometric band while 
dithering a few arcmins. The total exposure time 
was 1800 s (300 s x 6) per band. 

2.2. Data processing 

The data reduction was performed in a stan- 
dard manner, including overscan subtraction, bias 
subtraction, and flat-fielding with dome-flat, us- 
ing the IRAF software packag^]. The flux cali- 
bration was obtained from nearly 1,000 stars in a 
range of 15 - 19 mag which are listed in the Sloan 
Digital Sky Surve y (SDSS^I) Data R elease 8 pho- 
tometric catalog (jAihara et al.ll201ll) . The SDSS 



magnitudes (psfMag) are trans ferred to the BVR 
magnitudes (| Jester et al.ll2005lh and then corre- 
lated against instrumental magnitudes to derive 
the zero-points. The SEx tractor software package 
(jBertin fc Arnoutsl fl996l) was used with a circu- 
lar aperture of 8" radius. The uncertainty in the 
zero-points is estimated to be 0.01 - 0.02 mag. The 
systematic errors introduced through the conver- 
sion of magnitudes are approximately 0.03 mag. 

2.2.1. Masking of stars and galaxies 

In order to obtain the diffuse component, we 
masked stars and galaxies in the images. Our 
masking procedure is as follows. First, we de- 
tected all the discrete sources that exceed 1.5 a of 
the local background in more than three contigu- 
ous pixels using SExtractor with parameters DE- 



1 Kiso observatory is operated by the Institute of Astronomy 
of the University of Tokyo 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 

3 Funding for SDSS-III has been provided by the Alfred 
P. Sloan Foundation, the Participating Institutions, the 
National Science Foundation, and the U.S. Department 
of Energy Office of Science. The SDSS-III web site is 
http: / /www.sdss3.org/ 
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Fig. 1. — Optical field superimposed on the 
IRAS/DTRBE 100/im map. The solid square rep- 
resents the observed field. The brightest point 
of the lOO^m map corresponds to the center 
of CO emission from molecular cloud MBM32 



(Ma gnani et alJ ll985). 



TECT_ THRESH — 1.5 and DETECT _M1N AREA 
= 3. Second, from the detected bright stars, we 
derived point spread functions (PSFs) by fitting 
a Moffat function I = J [l + (r/a)- 2 ]^ (iMoffat 
1969), where Iq is the intensity at the center of the 
stellar image and a and /3 are the fitting parame- 
ters depending on the seeing. Figure [2] shows the 
radial profiles of bright unsaturated stars along 
with the best-fit curves. Adopting these best-fit 
curves as the PSF templates, every unsaturated 
star is masked with a circle with the radius where 
the count level is equal to 1/10 cr of the local back- 
ground count level. Saturated stars, magnitudes 
of which are supposed to be underestimated, are 
masked with a sufficiently large radius such as five 
or ten times larger than the masking radius es- 
timated by using the flux value and the PSF for 
an unsaturated star. The unmasked and masked 
images are presented in Figure [3] and Figure 01 re- 
spectively. The fractions of masked areas in B, g, 
V, and R images are 17%, 31%, 19%, and 35%, 
respectively. 

2.2.2. Aperture correction 

In the flux calibration of the surface bright- 
ne ss of extended sources, as e xtensively disc ussed 
bv lBernstein et al.l (|2002T l and iMattilal (|2003h . one 
must compensate for the flux that is scattered and 
diffracted into the aperture. Here we introduce 
the aperture correction factor, Ta, which is the 
fraction of the flux from a point source within the 
aperture. The flux from the point source, F(X), 
in units of Jy is given by 



(i) 



where C(A) is the signal counts within the aper- 
ture in instrumental units (ADU), Q(X) the sensi- 
tivity function in units of Jy ADU -1 . Because 
F(X) is known and C(A) is measured, we can 
only derive Q'(X)= Q(X) /T A (A), but not Q(X) 
and Ta (A) independently. In case of an extended 
source, no flux is lost; the flux lost from the aper- 
ture 1 — Ta(X) is compensated for by the flux scat- 
tered in, from the outside the aperture. Therefore, 
the brightness of the extended source 5(A) in units 
of Jy sr _1 is simply given as 



5(A) = ^C(X) 



YWTa(x) 
n 



C(A), (2) 
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Fig. 2. — The scaled radial profiles of stellar im- 
ages. Dots represent the normalized fluxes of ap- 
proximately hundred stars for each band. The 
red lines are the best fit of a Moffat function; 
J = / [l + (r/«r 2 ]^. 



where is the solid angle of the aperture in units 
of steradians. In the above equation, 5(A) is ob- 
tained once Ta(X) is given. 

The flux loss calculated from the fitted Moffat 
funct i on is negligible; ?a(A) = 0.994. However, 
Kind (|197ll ) pointed out that outside the central 
part of the profile, which is represented by a Gaus- 
sian or Moffat function, there is a more slowly de- 
clining halo or aureole. The aureole has a profile 
of an inverse-square law that extends out from the 
central part by a factor of 1000 in angular dis- 
tance, and contains about 5% of the star's light. 
Similar result s are presented by several authors 
( Mattilafl"2003l and references therein) . Even higher 
aureole energ y fractions of ~ 10% are reported 
in some cases dCapaccioli fc de Vaucouleurs 1983t 



Uson et al.l 119911 ). Here, we adopt the aperture 
correction factor Ta — 0.95±0.05 as a compromise. 

3. Analysis and Results 

Figure [5] plots the intensity of diffuse opti- 
cal light, 5„(A), as a function of 100/im inten- 
sity, 5 !/ (100//m). 5 y (100/im) is taken from the 
i"iL45/DIRBE 100/im map where the zodiacal 
light, stars, and galaxies are subtracted. We re- 
duced the resolution of the optical images to the 
same resolution as the 7_fL45/DIRBE 100/im map 
(2.'372 x 2?372) by averaging approximately 9,000 
optical pixels. Considering the seeing induced cor- 
relation, approximately 10 optical pixels consti- 
tute one independent data point. Thus, in the 
course of smoothing, the optical measurement er- 
rors are reduced by a factor of 30 as if averaging 
900 independent data points. The mean one a 
error per smoothed pixel, <r(A), of the smoothed 



Table 1: Observations and errors. 





^center 


a(A)t 


Obs. date 


Band 


W 


(kJy sr" 1 ) 


CUT) 


B 


0.44 


0.48 


2011-02-04 


9 


0.49 


0.39 


2011-04-03 


V 


0.55 


0.79 


2011-02-21 


R 


0.65 


0.39 


2012-02-21 


IRAS/BIRBE 


100 


0.3 xlO 3 





^Typical standard deviation for a smoothed 2f372 X 2f 372 pixel. 
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Fig. 3. — The optical images in 5-band (a), g-band (b), F-band (c), and _R-band (d). North is up, east to 
the left. The held is approximately 45' x 40', slightly smaller the field of view of instrument due to using 
dithering. 
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Fig. 4. — Same as for Figure 3, but for the masked optical images. 
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star-masked images is summarized in Table [TJ Fi- 
nally, we obtained 385, 336, 384, and 309 indepen- 
dent data points at B, g, V, and R, respectively. 

Figure [5] shows that S V (X) clearly correlates 
with S u (l00fj,m) at all the bands. We now fit to 
the data a linear function defined as: 



S„(A) = o(A) + 6(A)S*„(100^m) 



(3) 



where the slope parameter at wavelength A is 6(A) 
= AS V (X)/ AS u (100fim), and the constant param- 
eter a(A) represents components independent of 
the Galactic diffuse optical light (i.e., atmospheric 
airglow, zodiacal light, and any other light such 
as the COB). Note that A is B,g, V, or R for our 
broadband system, and a(A),6(A) > 0. 

To perform a minimum \ 2 analysis with the ef- 
fective variance method ( Orear 1982f) which can 
take into account both x and y errors, we need to 
know the standard deviation, er(100/im), which is 
not explicitly estimated by SFD98. We thus car- 
ried out simulations by generating noiseless data 
points which follow equation ([3]) with given a(A) 
and 6(A), and adding Gaussian random errors to 
these points. The noiseless 100/im data range from 
MIN(100/zm) to 8 MJy sr" 1 with a frequency dis- 
tribution oc ^(lOO/im) - ". At the low intensity 
end of 5j/(A) in Figure the correlation appears 
to be truncated suddenly because of the large scat- 
ter of the 100/Ltm data. The center of the truncated 
boundary along the S v (100fim) axis is located 
around 1.8 MJy sr -1 . So, we set MIN(100/zm) = 
1.8 MJy sr -1 . This was confirmed and justified by 
changing MIN(100/im) in the simulations. While 
cr(A) is taken from Table [IJ er(100^m) is a param- 
eter to be determined in the simulations. We find 
that a combination of n = 2 and tr(100/^m) = 0.3 
- 0.4 MJy sr" 1 reasonably reproduces the correla- 
tion and the frequency distribution as a function 
of ^(100/im). 

The bottom panel of Figure [5] shows the sim- 
ulated data at the g-band along with the linear 
function recovered from the y 2 minimum analysis 
as well as the assumed linear function to generate 
the mock data points. The top and middle panels 
plot the observed data along with the recovered 
linear function. a(A) and 6(A) are listed in Table 

m 



4. Discussion 

4.1. Variations in slope parameter 6(A) 

To explore the variations in slope parame- 
ter 6(A), we compiled 6(A) values from the lit- 
erature into a sample which is given in Table 
[H Figure [6] compares the data in this sample 
with model predictions. For the data points in 
the sample, the IRAS / DIRB E 100 /jm maps 



are u se d in Witt et al. ( 20081) . iMatsuoka et al 



(|201ll ). iBrandt fc Drainel (I2012I). and th is work , 
while early work, namely Laureiis et al. (1987). 
Guhathakurta fc Tysonl (|l989l ). |Palev et all (|l99lT ) 



and IZagurv et al. I (|l999h used the original IRAS 
100/im maps. The optical data are a collection 
of heterogeneous samples which were taken with 
different techniques and analyzed using differ- 
ent methods. The broadband data taken from 
the literature are: fields towards diffuse cloud 



Lynds 1642 studied bv lLaureiis et al.1 (|19871 ). four 
diffuse clouds (|6| > 2 5 °) at Bj, R, and / by 



iGuhathakurta & Tvsonl (jl989h . a co mpact high 
latitude source at B, V, and R by IPaley etall 



1991 ), a molecula r cloud MCLD123. 5+24.9 by 



Zagurv et al.1 ( 19991 ). and four high- latitude molec- 



ular clouds (MBM 30, 32 41A , and 41D) at B, G, 
R, and I bv I Witt et all (l2008h . The data studied 
by IMatsuoka et al.l ( 2011 ) were taken beyond the 
zodiacal dust zone with the Pioneer 10/11 Imag- 
ing Photopolarimctcr(IPP) at B and R in regions 
with S v (100/im) < 3 MJy sr -1 at Galactic latitude 
|6| > 35°, covering abo ut a quarter of th e whol e 
sky. The spectrum by IBrandt fc Draine ( 2012 ) , 
labeled "Full Sky Continuum" in their Figure 3, 
are obtained by analyzing the optical spectra of 
nearly 90,000 blank-sky spectra from the Sloan 
Digital Sky Survey (SDSS), the weighted regional 
coverage of which is concentrated at intermediate 
Galactic latitudes |6| = 30° ~ 45°. Their spec- 
trum has a large uncertainty in the flux calibra- 
tion due to their analyzing method. The spectra 
plotted in Figure [6] are scaled with their preferred 
bias factor of 2.1. In the case where more than one 
cloud is published in a single paper, the average 
6(A) value is plotted with the standard deviation. 
As can be seen, 6(A) varies by a factor of 3 - 4. 

In the following subsections, we explore possi- 
ble causes of the large 6(A) scatter by examining 
the effects of optical depth, dust albedo, dust tem- 
perature and the forward-scattering characteristic 



8 





315 




i i i 




simulation 




310 












305 


. ', -si".."! i * 

* ** 'jr^Xvt • 








HH ■ 


300 


■ ■ ■ ■ 





1 2 3 4 5 6 7 
S„(100^m) (MJy/sr) 



Fig. 5. — Correlation of the intensity of the diffuse optical light against 100/im intensity. The top and 
middle panels show the measured data. The red solid lines represent the linear functions recovered from the 
X 2 minimum analysis, while the red dashed lines represent the quadratic functions. In case of a quadratic 
function, the fit is limited to the data points with S v (100fim) > 2.0 MJy sr _1 (right to the vertical dash 
line); all the data points are used to fit the linear functions. The bottom panel shows the mock data. The 
black dashed line represents the linear function assumed in the simulations, while the red line represents the 
recovered function. g 



Table 2: DGL model parameters. 



S„(A) = a(A) + b(A)ff„(100/im) S„(A) = a Q {\) + b LJ (X)S v (100 f im) + c Q {\)S„(W0um)' 2 



Band 


a(A) 
(kJy sr" 1 ) 


6(A) 
(xlO -3 ) 


a«(A) 
(kJy sr" 1 ) 


6* (A) 
(xlO -3 ) 


c«(A) 
([kJy sr" 1 ]" 1 x 10- 6 ) 


B 


195.23 ± 0.12 


1.61 ± 0.11 


194.28 ± 0.50 


2.17 ± 0.28 


-0.08 ± 0.04 


9 


299.16 ± 0.16 


2.25 ± 0.14 


298.00 ± 0.60 


2.80 ± 0.34 


-0.07 ± 0.04 


V 


585.27 ± 0.46 


4.00 ± 0.28 


583.67 ± 1.02 


4.74 ± 0.56 


-0.11 ± 0.07 


R 


467.06 ± 0.22 


3.37 ± 0.21 


463.61 ± 0.91 


5.23 ± 0.47 


-0.24 ± 0.06 



of dust grains coupled to the non-isotropic inter- 
stellar radiation field (ISRF). 

4- 1.1. Optical depth and dust albedo 

The correlation in Figure [5] may appear to be 
linear. However, as we discuss below, linear cor- 
relations would not be expected in the case where 
the optically thin limit is not applicable, and sig- 
nificant changes in the optical depth along the 
sightline would cause large variations in 6(A) from 
cloud to cloud. 

Figure[7]depicts 6(A) as a function of S™ e (100/mi) 
in two ranges of the effective broadband wave- 
lengths, namely the blue range(0.43 - 0.48 /xm) 
and the red range (0.62 - 0.66 /mi). S^ 6 (100/im) 
refers to the representative average of the low- 
est and highest values of the (10G/im) range 
in which 6(A) are derived by assuming a lin- 
ear function. For the three c louds published in 
Guhathakurta & Tvsonl (|l989h . the plotted 6(A) 
values refer to individual clouds, not averaged 
values. The figure clearly shows that the 6(A) de- 
creases as S" ve (100fixn) increases, indicating that 
a significant portion of the reason for the large 
6(A) scatter seen in Figure [6] can be attributed to 
this 6(A) - S™ e (100/im) anti-correlation. 

Accordingly we propose a correction to equa- 
tion © as 6(A) = 6°(A) + 6 1 (A)S , „(100/im), or 
S„(A) = a(A) + [6 (A)+6 1 (A)S*„(100Aim)]S' l/ (100/xm). 
Then, S V {X) has a quadratic term as follows: 

S v {\) = a Q (A)+6 Q (A)A(100/im)-|-c Q (A)&(100/im) 2 , 

(4) 

where b Q {\) > and c Q (X) < 0. While the cor- 
relation in Figure [5] appears to be linear, the ac- 
tual fact is that the correlation is non-linear, but 
this is hidden by the prevailing noise. Since the 
quadratic term c < ^(A)S' l/ (100/j,m) 2 is negative, this 
would explain why the correlation of the Pioneer 
data deviates below the best-fit linear function at 



the b right end of 100/tm emission (jMatsuoka et al 
20111 ). We performed a minimum x 2 analysis with 
the effective variance method to fit the quadratic 
function to our data. As can be seen in Table [2] 
and Figure [SJ the derived constant parameter of 
a quadratic function a®(X) is smaller than that 
of a linear function a(A) at all the photometric 
bands. Thus, this negative quadratic term would 
be important for the separation of the COB from 
the diffuse optical light, reducing the DGL inten- 
sity at which 5V (100/im) is equal to the Cosmic 
100/im background intensity. 

The origin of the quadratic term could be at- 
tributed to the fact that our sample directions are 
not optically thin, but translucent or thick. Fig- 
ure M illustrates the relation between reddening 
E(B — V) and 100/im emission at the peak of 
100/tm emission for our clouds, where E(B — V) is 
taken from the IRAS /DIRBE reddening map. If 
typical extinction curve (Ry = Ay/E(B — V) = 
3.1) is assumed, the 100/im brightness per unit of 
visual extinction is S v (100^im)/Av = 8-15 MJy 
gr -i ma g-i Our sample ranges from ^(lOO/im) 
= 3 to 18 MJy sr" 1 , corresponding Ay — 0.16 - 
2.0 mag. Thus, dust absorption along the sightline 
should be taken into account. 

To illustrate the effect of extinction along the 
sightline, as fully discussed in the Appendix, we 
consider a simple model where a plane-parallel 
dust slab is illuminated by the monochromatic 
starlight from the back. The 6(A) is expressed as 



6(A) oc exp[— (1— u;)t][1— exp(— u;t)]/{1— exp[— (1— cj)t]}, 

(5) 

where u is the albedo. 

Equation implies that 6(A) depends on the 
optical depth and the albedo, but not on the 
ISRF intensity. In Figure [71 we plot equation ([5]) 
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Table 3: Correlation slopes 6(A) for high-latitude clouds or fields. 



Reference 
Clouds or fields 


b(\y 

(xlO -3 ) 


^center 

(pm) 


5^(100/im) range* 
(MJy sr" 1 ) 


b(R)/b(B) 


This work 










MBM32(/ = 147,6 = 41) 


1.6 ± 0.1 

2.3 ± 0.1 
4.0 ± 0.3 

3.4 ± 0.2 


0.44 
0.49 
0.55 
0.65 


1 - 7 


2.1 


Laureijs et al. 1987 










ijynasro4z^t — ^rr, o — — oo ) 


z.z it u.o 


n at 


1 1 Q 




Guhathakurta et al. 1989 










irl(Z = 174, 6 = -42) 


0.4 
1.2 


0.45 
0.65 


7 - 18 


3.2 


iv2(l = 235, 6 = 37) 


1.1 
2.2 


0.45 
0.65 


3 - 7 


2.0 


ir3(Z = 38, b = 45) 


2.6 

A A 

6.0 


0.45 
u.uo 
0.90 


3 - 9 


1.7 


Paley et al. 1991 










CTI cirrus cloud(Z = 104, 6 = -32) 


4.5 ± 0.3 
8.5 ± 0.8 
1 1 n 4- i n 


0.44 
0.55 

u. / U 


2.5 


2.3 


Zagury et al. 1999 










MCLD123.5+24.9(Z = 124, 6 = 25) 


2.4 ± 0.1 

A Q _l_ n i 

3.0 ± 0.1 


0.45 
u.uo 
0.90 


5 - 8 


1.8 


Witt et al. 2008 










Derived from 4 clouds* (Qft < / < 1 47 38 < h < 41 1 


2.1 ± 0.1 
3.3 ± 0.4 

4.2 ± 0.7 
2.7 ± 0.7 


0.46 
0.53 
0.63 
0.83 


1-8 


2.0 


Matsuoka et al. 2011 










1/4 of the entire sky (|6| > 35°) 


2.1 ± 0.1 
4.6 ± 0.1 


0.44 
0.64 


1 - 3 


2.2 


Brandt et al. 2012 










SDSS spectra (30 < b < 45)** 


1.2 ± 0.1 

2.3 ± 0.1 


0.45 
0.65 


1 - 10 


1.9 



^The correlation slope is defined as a linear function of S V (X) = a(\) + 6(A)5, y (100/^m) 

range refers to the range in which fe(A) is derived assuming a linear function. 
*MBM30(Z = 142.2, b = 38.2), MBM32(Z = 147.2, b = 40.7), MBM41A(Z = 90.0, b = 39.0), MBM41D(Z = 92.3, b = 37.5) 
"The weighted regional coverage of the analyzed nearly 90,000 blank-sky spectra is concentrated at intermediate Galactic latitudes 
\b\ = 30° ~ 45°. 
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this work (MBM32, S„(10qam) < 7 MJy sr" ) 
Laureijs + 1987 (Lynds1642, S„(10C|«m) < 13 MJy sr ) 
Guhathakurta+1989 (4 clouds, S„(100^m) < 18 MJy sr" 1 ) 
Paley+1991 (1 cloud, S v (10q«m) < 2.5 MJy sr" 1 ) 
Zagury+1999 (MCLD1 23.5+24.9, S v (100^m) < 8 MJy sr 1 ) 
Witt+2008 (MBM30, 32, 41 A, 41 D, S„(10q«m) < 8 MJy sr J 
Matsuoka+2011 (|b| > 35degree, S„(100^m) < 3 MJy sr" ) 
Brandt+2012 (SDSS spectra, S^IOO/^m) < 10 MJy sr" ) 
Brandt+2012 (model ZDA04 + de-reddened MMP83) 
Brandt+2012 (model WD01 + de-reddened MMP83) 



_i i i i i i i i i_ 



_1 I I I I I I I L 
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.1 I I I I I I I L 



j i i i i i i i i_ 
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0.6 0.7 
Wavelength (fim) 



0.8 



0.9 



Fig. 6. — The correlation slopes 6(A) = AS , ;y (A)/AS' I y(100^m) as a function of wavelengths. 
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with a single constant value of proportionality^!, 
X, both for the left and right panels, assuming 
S u (lQ0(im)/A v = 10 MJy sr" 1 mag" 1 . X is set 
to an appropriate value by eye. The three lines 
depict the relation for albedo ujy of 0.6, 0.7, and 
0.8 from bottom to top. Figure [7] suggests that a 
range of Ay = 0.16 - 2.0 mag could change the 
6(A) strength by a factor of 2 - 3 while keeping ui 
constant. It also indicates that the Ay range alone 
cannot explain the full range of the 6(A) scatter. 
If the albedo becomes greater as the optical depth 
becomes smaller, then the full range of the 6(A) 
scatter could be explained in terms of the opti- 
cal depth effect coupled to variations in albedo. 
However, this case is not likely, because small-size 
dust grain s, which have smaller albedo than large 
ones (e.g. |Prainell2003|) . are expected in smaller 
optical depth regions according to the theory of 
grain growth. 

Two smooth curves in Figure [HI are synthesize d 
spectra for ry = 0.15 bv lBrandt fc Drainel (120121) 



adopt ing the ISRF spectrum by iMathis et al 
() 19831 ) combined with a simple pl ane-parallel 



galaxy mod el and two dust models from | Zubko et al 
(|2004h and IWeingartner h Drainel (|200l[ ) fhere- 
after ZDA04 and WD01 models). The models 
underestimate 6(A) by a factor of 2. There would 
be two possible explanations for this. One is a de- 
ficiency in UV photons compared to the assumed 
spectrum of the ISRF; the dust temperature is low 
and 6(A) thus increases where the dust grains are 
illuminated by a sm all number of UV photo ns (see 
Figures 9 and 10 in lBrandt fc Drainell20l3) . The 
other is that the albedo is greater than assumed in 
the models, because Figure [7] and equation ([5]) im- 
ply that the 6(A) strength is doubled if lu increases 
by 20%; for example, uy to be 0.8 inste ad of a 
stand ard Milky Way value of 0.67 at V (jPrainel 
l2003h . If this is the case, more large grains are 
needed than in the ZDA04 and WD01 models, 
because large grains are efficient scatters of light. 

4-. 1.2. Dust Temperature 

Another possible cause of the large 6(A) scat- 
ter is variations in dust temperature. Higher 
dust temperature increases the intensity of the 
100/im emission and, as a consequence, it would 
decrease the 6(A). The temperature dependence 




5 10 15 

S„ (lOQam) (MJy/sr) 



20 



Fig. 8 — E(B -V) as & function of S„(100/im) for 
the peak of 100/im emission in the cloud. The legend 
for the points is the same as in Figure [6] 



of the resulting 100/im r adiation is ~ T 5 - 8 at 18K 
(|Brandt k. Drainell2012l ). However the real situa- 
tion is not so simple. The temperature of 100/im- 
emitting dust grains, i.e. large grains, is deter- 
mined by the balance between incoming energy 
from the illuminating radiation to the grains and 
the outgoing infrared thermal energy. One can 
increase the temperature of the grains by increas- 
ing the UV component of the ISRF while leav- 
ing the optical intensity of the ISRF the same. 
That would leave the optical DGL intensity the 
same but increase the 100/im intensity. However, 
if the optical intensity of the ISRF changes, it also 
changes the intensity of the DGL. Thus, the incre- 
ment of 6(A), e.g the ratio of increment of the DGL 
and that of the 100/im emission, depends intri- 
cately on the spectral ener gy distribution and the 



4 b(A) =Xexp[-(l-a))r][l-exp(-a)T)]/{l-exp[-(l-£j)r]} 



non-isotropy of the ISRF. iLehtinen et al.1 (|2007) 
showed that the dust emissivity, which is measured 
from the ratio of far-IR optical depth to visual 
extinction r (far-IR) /Ay, decreases as a function 
of dust temperature, and the higher temperature 
dust does not necessarily mean the higher surface 
brightness at 100/im. They concluded that the 
emissivity variations are caused by the variations 
in the absorption cross section of dust in the far- 
IR, or the grain size distribution. 

In our sample, we do not see any convinc- 
ing correlation between 6(A) and dust tempera- 
ture T as shown in Figure [9] The plotted dust 
temperatures are taken from the IRAS/DIRBE 
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Fig. 7. — The left panel plots the correlation slopes 6(A) as a function of S , °" e (100/xm) for the blue range (0.43 to 
0.48 frni). S , "" e (100/xm) refers to the representative average of the lowest and highest values of the ^(lOO^im) range 
in which 6(A) is derived by fitting a linear function. Spectroscopic data are averaged over this range. The right panel 
is the same as for the left panel but for the red range (0.62 to 0.66 um). The le gend for the points is the same as in 
Figure [SI except for the crosses which are the averaged spectroscopy data from lBrandt Draind (|2012f ). The lines 
represent the relation expected from equation (JSj for albedos uiv = 0.6, 0.7, and 0.8 from lower to upper. 



temperature map. The T 5 8 dependence of the 
100/im intensity causes an increase by a factor 
of about 1.7 over a temperature range from 16K 
to 17. 5K shown in the data of Figure |9] How- 
ever, the ratio of r(far-IR)/Av declines with in- 
creasing temperature, and the effective increase in 
S^lOO/im) would become lower than 1.7. This 
may explain why the data of Figure [9] do not 
lead to any discernible trends. As a caveat, it 
should be noted that the IRAS /DIRBE temper- 
ature map is severely limited by the angular res- 
olution of DIRBE of ~ 1 degree. Thus, this map 
is not capable of correctly reflecting the tempera- 
ture variations across cloud structures with scales 
much smaller than the DIRBE resolution. Further 
exploration of the 6(A) dependence on the dust 
temperature requires far-IR data with higher res- 
olution such as those from the Akari all sky survey. 
Such analysis is beyond the scope of this paper. 

4-1-3. Forward-scattering characteristic of dust 
grains 

Another possible cause of the large 6(A) scatter 
is the effect of the forward-scattering characteristic 
of dust grains. The ISRF is non-isotropic because 
of the stellar concentration toward the Galactic 
disk. The scattering efficiency in the direction to- 
wards the observer is coupled to this non-isotropy 



through the forward-scattering phase function of 
scattering grains. This effect is maxi mized for the 
geo metric configuration explored by lJural (Il979l) 
and IStarkl f)l995h . The latitude dependence of 
surface brightness of a cl oud at high latitude 
was st udied, employing the iHenvev fc Greenstein 
([19411 ) approximation for the scattering phase 
function. In this model, it is assumed that a 
cloud is illuminated by an infinite homogeneous 
disk. The observer is located in the plane, and 
the cloud is lo cated above the plane. The numer- 
ical results by lJura (Il979l) show that the surface 
brightness of a cloud S is written as 



S cx 1 — l.lg \J sin |6| 



(6) 



where g is the asymmetry factor of the phase func- 
tion and 6 the Galactic latitude. For strong for- 
ward scattering, g ~ 1; for isotropic scattering, 
g ~ 0. This equation shows that the brightness of 
the scattered light is significantly changed with the 
latitude of the cloud. Adopting a typical asym- 
metry factor g ~ 0.75, the intensity of the scat- 
tered light at 6 = 30° is 2-3 times larger than 
that of at 6 = 90°. This model assumes that all 
radiation sources, i.e. stars, are located just in 
a plane, however, this is not the case. In fact, 
the mean v ertical height of high- latitude clouds is 
~ 150 pc ( Magnani et al. 19961) and that of the 
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Fig. 9. — The left panel plots the correlation slopes 6(A) as a function of T<i ust for the blue range of the effective 
broadband wavelengths. The right panel is the same as for the top panel but for the red range. The legend for the 
points is the same as in Figure [6] The horizontal error bars are represent their temperature ranges among the clouds. 



stellar disk is 



300 p c (|Bahcall fc Soneiralll98d: 
Gilmore fc Reidl 19831). Thu s, the geometry used 



m the model of Jural ( 19791) is extreme. It would 
be more natural to assume that a significant por- 
tion of stars are located behind the high-latitude 
clouds. 

To examine the effectiveness of the forward- 
scattering characteristic of dust grains, we present 
Figure [TUl illustrating 6(A) as a function of Galac- 
tic latitude 6| (upper panels) and longitude \l\ 
(lower panels). In upper panels, we may see a 
correlation between the 6(A) and the Galactic lat- 
itude |6 in the blue range, but no correlation in the 
red range. There are no 6(A) dependency on the 
longitude \l\ at all. For comparison, equation © 
with g — 0.75 is plotted in upper panels of Figure 
HUl The absolute value of the equation is deter- 
mined by eye to fit the observed data. The 6(A) 
variation in the blue range seems to be partly ex- 
plained in terms of the forward-scattering. Then, 
why are there no correlations in the red range? A 
possible cause could be a difference in the distri- 
bution of blue stars (i.e., effective emitter of blue 
light) and red stars (i.e., effective red light emit- 
ter). Blue young stars are more concentrated to- 
wards the Galactic plane than red old stars and 
the g value at the blue rang e is greater t han the 
red range in dust models (e.g. lDrainell2003l ). Thus, 
the effect of the forward-directed scattering phase 
function would be more efficient in the blue range 



than in the red range. 



4.2. Extended red emission 



Extended red emission (ERE) is observed as 
an excess m the 5000-9000A spectral range in the 
ISM. It is suggested that ERE is the result of 
an inter action of far-UV photons with interstel- 
lar dust (IWitt fc Schildlll985UDarbon et al.lll999t 



Witt et al l 120061) . ERE is observed in reflection 



nebulae Jwitt et all Il983 IWitt fc Schiidl fl98 



Witt fc Borosonll9 90), HII regions (Perrin & Siva: 



1990 



1992), and planetary nebulae ( Furton fc Witd 
19921 ). In the diffuse ISM, such as high 



latitude clouds, some authors favor the pres- 
ence of ERE in the diffuse optical lightn (e.g. 



i Guhathakurta fc Tvson1ll989l: iGordon et aL| [l9M 
IWitt et al.ll2008HMatsuoka et al.ll201ll). while oth 



ers suggest the ab sence of ERE (jZagurv et al 



19991: IZagurvl 120061) based on broadband observa 



tion s. The c l oud ob served by IZagurv et al. I (ll999l) 
and Zagurvl (|2006l ) is illuminated by the nearby 
star Polaris (F7 lb). The optical surface bright- 
ness of the cloud is dominated by scattered light 
of Polaris, but this star is not blue enough to emit 
sufficient far-UV photons for ERE excitation (pri- 
vate communications with the anonymous referee; 



5 In this section, we use the term "DGL" as the scattered 
component and "diffuse optical light" as the total optical 
component including the DGL and any other emissions. 
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Fig. 10. — The upper panels plot the correlation slopes 6(A) as a function of Galactic latitude |6|, and the lower 
panels the correlation slopes b(\) as a function of Galactic longitude \l\. The left and right panels are for the blue 
range of the effective broadband wavelengths and for the red range, respectively. The line solid curve represents the 
prediction from equation ||B). which is derived by assuming models where the effectiveness of the forward-scattering 
characteristic of dust grains is maximized. The legend for the points is the same as in Figure [6] The horizontal error 
bars are represent coordinate ranges. 
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see section 6.4 of iLehtinen fc Mattilal (|2013[) for 
further discussion). Thus, it would be natural 
for them to miss ERE. Carrying out spectroscopy 
of the diffuse optical light toward a high latitude 
cloud, Szomoru fe Guhathakurtal dl998h claimed 
to detect ERE peaking at 6000 A, which is blue- 
shifted compared with 8000A for HII regions and 
7000A for reflection nebulae and planetary nebu- 
lae. On the othe r hand, the more r ecent spectro- 
scopic study by iBrandt fc Draind ( 20121 ) cannot 
find any evidence of ERE. 

There is some confusion in the DGL models, 
which are sometimes used for the claim of ERE 
detection; ERE is detected, if the observed diffuse 
optical light is redder than predicted by the model, 
and vice verse. Here we introduce the b(R)/b(B) 
value, which is a measure of the redness of the 
diffuse optical light. As indicated in Tabled mea- 
su red b(R)/b(B) values a re typically ~ 2. Models 
bvlWitt fc Schildl (ll985l).lGuhathakurta fc Tvsonl 
(|l989l) . and lBernstein et al.l (|2002h predict b(R)/b(B) 
1.4, which is bluer than the measured values 
b(R)/b(B) ~ 2, supporting the presence of 
ERE in the diffuse optical light. In contrast, 
iGordon et ail (jl998l) predicted b{R)/b(B) ~ 1.9, 
using the Witt-Petersohn DGL model(WP model; 
Witt fe Petersohnlll994l: IWitt et al.lll997t ). an ad- 



vanced Monte Carlo multiple scattering model uti- 
lizing the actual Galactic radiat ion field and real- 
i stic d ust scattering properties. IBrandt fc Draine 
(|2012l ) also predicted b{R)/b(B) ~ 1.9 - 2.0 
using their scattering DGL modcls(ZDA04 and 
WD01). These are consistent with the observed 
values without ERE. However, the observed val- 
ues deduced from PIONE ER10/11 dat a from 
the high-latitude regions by IGordon et al. (Il998l) 
are b(R)/b(B) ~ 3.2, which is higher than other 
observed values, and still redder than the models. 

It is emphasized that 6(A) in this work on cloud 
MBM32 in Figure [B] has an excess at the U-band 
and favors ERE in the diffuse optical light; b(X) 
rises from B to V faster than the models, seems to 
peak around 6000 A and decreases towards long 
wavelengths. Such a characteristic is expected 
from the models in which scattered DGL is com- 
bined with ERE. 



5. Summary 

We have conducted B, g, V, and i?-band imag- 
ing in a 45' x 40' field containing part of the high 
Galactic latitude translucent cloud MBM32, and 
correlated the intensity of the diffuse optical light 
S U (X) with that of 100/zm emission S^lOO^m). 
5^(100/im) ranges from 1.4 to 7.0 MJy sr _1 in the 
field. A minimum \ 2 analysis is applied to fit a 
linear function to the measured correlation and to 
derive the slope parameter 6(A) of the best-fit lin- 
ear function. 

Compiling a sample by combining our 6(A) and 
published ones, we showed that the 6(A) strength 
varies from cloud to cloud by a factor of 4. We 
also found that 6(A) decreases as ^(lOO^m) in- 
creases in the sample. To incorporate this 6(A) - 
5^(100/im) relation, we proposed that a non-linear 
correlation including a negative quadratic term of 
S^lOO^m) 2 be fitted to the measured correlation. 

We explored possible causes of the large 6(A) 
scatter by examining the effects of the optical 
depth, dust albedo, dust temperature and the 
forward-scattering characteristic of dust grains 
coupled to the non-isotropic interstellar radiation 
ficld(ISRF). The variation of optical depth, which 
is Ay = 0.16 - 2.0 in the sample, can change 6(A) 
by a factor of 2 - 3 while keeping the albedo, ui, 
constant. This indicates that the Ay range alone 
cannot explain the full range of the 6(A) scat- 
ter. Our finding of a correlation between the 6(A) 
and Galactic latitude |6| in the blue range (0.43 
- 0.48/Ltm) would imply that there would be some 
contribution to the 6(A) scatter from the effect of 
the forward-scattering characteristic of dust grains 
coupled to the non-isotropic interstellar radiation 
field(ISRF). 



The models by IBrandt fc Draind (|2012j ) under- 
estimate 6(A) by a factor of 2. There would be two 
possible explanations for this. One is a deficiency 
in UV photons compared to the assumed spectrum 
of the ISRF, because fewer UV photons result in 
lower dust temperature, which increases 6(A). The 
other is that the albedo is greater than assumed 
in the models. If this is the case, a greater number 
of larger grains are needed in the models, because 
large grains are more efficient scatterers of light. 

Our 6(A) spectrum favors the presence of ex- 
tended red emission(ERE) in the diffuse optical 
light; 6(A) rises from B to V faster than the mod- 
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els, seems to peak around 6000 A and decreases 
towards longer wavelengths. Such characteristic 
is expected from the models in which the DGL is 
combined with ERE. 
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A. Radiative transfer in a dusty slab 

Here we consider a simple model where a plane-parallel dust slab is illuminated from the back by monochro- 
matic starlight at optical wavelength A. We measure optical depth r inward toward the plane-parallel slab. 
Light is produced along a sightline by the scattering of the illuminating starlight and from scattered light in 
the slab. Here, we assume complete forward-scattering for simplicity. The contribution of an element with 
optical depth from t to t + dr to the scattered light and the starlight are given by 

dlu,aca(^) = [Iu,star(r) + I u , SC a{T)]ujdT - I„, S ca{T)dT (Al) 
dI V: star(T) = -Iu,star(T)dT (A2) 

where I VtS ca( T ) an d Iv,star{T~) are the specific intensity of the scattered light and the starlight, respectively, 
U) is the albedo. Then the solution of I^sca(i~) is given by 

Iu,sca{r) = Iv.star (0)exp[-(l - u)t] [1 - exp(-wr)]. (A3) 

The far-infrared (IR) luminosity L(IR) comes from thermal emission from dust grains which absorb the 
starlight and the scattered light. The contribution of an element to L(IR) is given by 

dL v<IR (T) oc [Iv.starir) + I v<sca (t)] (1 - uj)dr. (A4) 

By integrating this equation, we have 

L v ,ir{t) oc l v , star (0){l - exp[-(l - u)t]}. (A5) 

The conversion of the total IR luminosity to the 100/mr luminosi ty is L(lOO^m) = 0.52 L(IR), which is 
applicable for models from 0.5 to 1.5 times the local ISRF intensity ( Brandt fc Draine l2012h . Thus we have 

6(A) oc exp[-(l - w)r][l - exp(-wr)]/{l - exp[-(l - w)r]} (A6) 

where 6(A) is the intensity of the scattered light at wavelength A relative to the intensity of 100/im emission. 
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